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Abstract. We have collected all of the published pho- 
tometry for GRB 990123 and GRB 990510, the first two 
gamma-ray bursts where breaks were seen in the light 
curves of their optical afterglows, and determined the 
shapes of their light curves and the break times. These 
parameters were used to investigate the physical mech- 
anisms responsible for the breaks and the nature of the 
ambient medium that the bursts occurred in. The light 
curve for GRB 990123 is best fit by a broken power law 
with a break 1.68 ± 0.19 days after the burst, a slope 
of ai = — 1.12 ± 0.08 before the break, and a slope of 
OL2 — —1.69 ±0.06 after the break. This is consistent with 
a collimatcd outflow with a fixed opening angle of #o ~ 5° . 
In this case the break in the light curve is due to the rel- 
ativistic fireball slowing to T « l/#o- The light curve for 
GRB 990510 is best fit by a continuous function with an 
early-time slope of a\ = —0.54 ± 0.14, a late-time slope 
of a.i — —1.98 ± 0.19, and a slow transition between the 
two regimes approximately one day after the burst. This 
is consistent with a collimated outflow with 9q « 5° that 
is initially radiative, but undergoes a sideways expansion 
that begins approximately one day after the burst. This 
sideways expansion is responsible for the slow break in the 
light curve. 

Key words: Gamma rays: bursts 



1. Introduction 

In the first few months of 1999 bright optical afterglows 
(OAs) were observed for two gamma-ray bursts (GRBs): 
GRB 990123 and GRB 990510. Extensive follow-up obser- 
vations were made at optical wavelengths for both GRBs, 
so it was possible to follow their light curves down to 
V w 28, with the deepest photometry being done with 
the Hubble Space Telescope (HST). These data, when put 



onto a common photometric system, can be used to con- 
strain the shape of each burst's light curve. The resulting 
shape parameters (the time of the break, and the slopes 
before and after the break) can be used to constrain physi- 
cal models of the OA, and the interaction of the relativistic 
ejecta from the initial blast with the surrounding medium. 

At 9:46:56.12 UT on 1999 January 23 the Burst 
and Transient Source Experiment (BATSE) detectors on 
the Compton Gamma Ray Observatory satellite detected 
GRB 990123 in the constellation Bootes. X-ray observa- 
tions from the Dutch-Italian satellite Beppo SAX localized 
the position of the burst to within 5' (Piro 199E ) and the 
Robotic Opti cal Transient Se arch Experiment (Akerlof & 
Akerlof et 



McKay 1999 



al. 



1999) 



detected an OA that 
= 8.86 ± 0.02 just 47 



reached a peak magnitude of V 
seconds after the BATSE trigger. Further optical observa- 
tions constrained the r edshift to betwee n 1.60 < z < 2.04 
(Andersen et al. 199E; Kulkarni et al. 1999a), and im- 



ages taken with the HST Space Telescope Imaging Spec- 
trograph (STIS) revealed a host galaxy with Vq — 24.25 
(Holland & Hjorth 1999). The burst had a peak fluence of 



(5.09±0.02) x 10" 4 erg cm~ 2 (Kippen 1999a), which places 



it in the top 0.3% of the BATSE fluence distribution. This 
corresponds to a total energy release, assuming isotropy 
and depending on the details of the adopted cosmology, 
of between « 3 x 10 54 and rs 4.5 x 10 54 erg, which is 
approximately the rest-mass energy of two neutron stars. 
This amount of energy is difficult to explain using merging 
compact objects since it would require that the entire rest 
mass of both objects be converted to gamma radiation 
with an efficiency of nearly 100%. 

GRB 990510 (the Anja burst) was detected by both the 
Compton Gamma Ray Observatory and BeppoSAX in the 
constellation Octans at 8:49:05.95 UT on 1999 May 10. 
Optical follow-up observations 3.5 hours after the burst 
revealed an OA with R m 17.5 (Axelrod et al. 
Vreeswijk et al. 
shift of z 



(1999) 



1999) 



found a lower limit for the red- 
1.619. A faint galaxy (V = 28.5 ± 0.5 ) has 
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been detected 07066 ± 07009 from the GRB (Bloom [2000 
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Fruchter et al 200C). The burst had a peak fluence of 
(2.56±0.09) x 10 5 erg cm" 2 (Kippcn |l999b| ), which places 
it in the top 9% of the BATSE fiuence distribution. This 
corresponds to a total isotropic energy release, depending 
on the details of the cosmology, of at least 10 53 erg, ap- 
proximately one order of magnitude less than the energy 
released in GRB 990123. 

The light curve for GRB 990123 became steeper ap- 
proximately two days after t he bur st (Castro-Tirado 
et al. 1999| , Kulkarni e t al. 1999a ), while Ha rrison 



et al. (|1999| ), Israel et al. ( |l999|) , and Stanek et al. ( |l999D 



found that the light curve of GRB 990510 became sig- 
nificantly steeper approximately 1.5 days after the burst. 
These breaks have been interpreted as evidence that the 
outflows from the bursts are collimated with opening an- 
gles of approximately 5°-10° (Sari et al. 1999, Castro- 
Tirado et al. 1999). Such models provide a solution to the 
so-called "energy crisis" of GRBs since, if GRBs are col- 
limated outflows, the total energy requirement drops by 
factor of between approximately 100 and 1000. This makes 
the energetics of GRBs consistent with the energetics of 
supernovae and merging compact objects. 

Before GRB 990123 GRB light curves were usually fit 
by power laws of the form /„(i) = kt a where k is a normal- 
ization constant and a is between approximately —1 and 
—2. The steepest temporal decay was a = —2.5 for th e OA 
associated with GRB 971227 (Djorgovski et al. |1998| ). The 
first evidence that the decay of the optical light curve from 
an OA varies with time came from GRB 990123, where a 
rapidly steepened from — 1.10 ± 0.03 to — 1.65 ± 0.06 with 
the break occurring 2.04±0.46 days after the burst (Kulka- 
rni et al. |!999aD - Har rison et al. fl!999| ), Israel et al. ( |1999| ), 
and Stanek et al. ( 1999] ) found that the optical light 
curve for GRB 9900510 exhibited a similar behavior, but 
with a slow transition between the two regimes. Harrison 



et al. (1999) used a four-parameter fitting function to find 
slopes of -0.82 ± 0.02 before the break and -2.18 ± 0.05 
after the break with the break occurring 1.20 ± 0.08 days 



after the burst. However, Stanek et al. (1999) found slopes 
of -0.76 ±0.01 before and -2.40 ±0.02 the break with the 
break occur ring 1.57±0.03 days after the burst; and Israel 
et al. ( |1999[ ) found -0.88 ± 0.03 before and -2.68 ± 0.13 
after with a break at 1.8 ± 0.2. The latter two groups 
used a different four-parameter fitting function from Har- 
rison et al. ( |1999D , and the fits formally disagree with each 
other. 

In this paper we collect the published photometry for 
the OAs associated with GRB 990123 and GRB 990510, 
examine three fitting functions for broken light curves, and 
use these fits to constrain physical models for the burst. 
We investigate the effects of using different fitting func- 
tions on determining the time of the break and the slope 
of the light curve before and after the break. Our goal is 
to determine reliable light curve parameters for each GRB 
and to use these parameters to distinguish between differ- 
ent physical models for the nature of each OA. Sect. 



discusses the three fitting functions that we used to deter- 
mine the break time, and the slopes of each light curve. 
Sect. |]presents our fits to the light curve of GRB 990123's 
OA and an estimate of the slope of the optical spectrum 
of the OA. Sect. ^ presents our fits to the light curve of 
GRB 990510's OA and an estimate of the slope of its op- 
tical spectrum. We discuss the effects of the choice of the 
fitting function on the determination of the break time 
and slopes of the light curves for each burst. In Sect. || we 
discuss various models for the decay of the light curves, 
and in Sect. ^ we derive some physical properties for the 
bursts based on our fits to the light curves. We have writ- 
ten these sections in such a way that they can be used 
with data from other GRBs. 

We have adopted a cosmology with a Hubble param- 
eter of Hq = 65 km s _1 Mpc -1 , and assumed a den- 
sity parameter of f2o = 0-2 and no cosmological constant 
(f^A = 0) throughout this paper. 

2. Fitting the Light Curves 

The choice of fitting function is arbitrary, but that choice 
can affect the determination of the slopes and break time, 
and thus can influence the interpretation of the light curve. 
For example, Jensen et al. ( 2000| ) find different break 
times and slopes for the decay of the optical light curve 
of GRB 000301C depending on whether a broken power 



law or the continuous function of Stanek et al. (1999) is 



fit to the data. We have chosen to fit the light curves for 
the OAs associated with GRB 990123 and GRB 990510 
with three fitting functions. The first of these is a broken 
power law, 



U(t) = 



U(tb)(t/t b ) ai , if t<t b 
fAt b )(t/t b ) a2 , if t>t b , 



(i) 



where f v (t) is the flux in /iJy t days after the burst, and 
tb is the time that the light curve made the transition 
from a slope of a\ to a slope of a-i. The slope goes as 
f u oc t ai before the break and /„ cx t a2 after the break 
with an instantaneous transition between the two slopes. 
This is the simplest form that can be fit to the data, and 
it can be reduced to two linear functions (i.e., In (f v {t)) = 
ai(t/tb) + \n(fv(tb)) before the break and In (/„(£)) = 
ct2{t/tb) ± In (/„(£;,)) after the break), which simplifies the 
fitting process. 

The second fitting function is, 



fu(h) 



where J 



{t/h 



(2) 



This is the continuous function 



used by Harrison et al. (1999). It is differentiable at t = 
tb, and is generally more flexible than Eq. (|l|) is. It is 
able to model the slow transition between the ct\ and oci 
regimes, which may be more physically realistic than the 
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instantaneous transition that is assumed by the broken 
power law. 

The third fitting function is, 



Table 1. The central wavelengths and photometric zero 



points for each band pass (from Fukugita et al. 1995). 



fu(t) = 



Vu{t b ){t/t b ) C 

i + (t/t b ) a2 -' 



(3) 



which is the same as that used by Stanek et al. ( [1999 ) 
and Israel et al. (1999), except that we have redefined 
and a2 to be consistent with our Eqs. ([!]) and (g). This 
function has the same advantages over Eq. (0) that Eq. (§) 
does. 

We wish to stress that none of these functional forms 
have been derived from any physical principles, and that 
they are not based on any model for the GRB or its OA. 
The three forms were selected solely because they allow us 
to parameterize the observed light curves for GRB 990123 
and GRB 990510. We used CERN's Minuit function min- 
imization package, and a chi-square minimization scheme, 
to simultaneously solve for the four free parameters (a±, 
ct2, tb, and f v (tb)) and to compute the formal la errors in 
the fit for each parameter. 

3. GRB 990123 Data 

3.1. The Photometry 

We collected the published optical observations of the 
OA associated with GRB 990123 from six sources: 
Castro-Tirado et al . ( p.999| ), Fruchter et al . fll999b[ ), 
Gala ma et al . ( |1999a| ), Holland & Hjorth Ql999| ), Kulkarni 
et al. ( |1999a[ ), and Fruchter et al. ( [1999c] ). These data con- 
sists of 88 observations in the Johnson B-, and V-bands, 
Kron-Cousins R-, and /-bands, and Gunn 7--band. The ob- 
servations were taken at 19 telescopes over a two month 
period after the burst. In cases where the same observa- 
tion was reported in two or more sources, and the quoted 
magnitudes disagreed with each other, we discarded both 
observations. This was done because we had no reason to 
favour one source over the other. In cases where a Gunn 
r-band magnitude was derived from an i?-band magnitude 
we discarded the derived Gunn r-band data and kept the 
original i?-band data. 

We converted the published magnitudes to fluxes 
using the photometric zero points, from Fukugita 
et al. (1995). These zero points, and the central wave- 



lengths for each band pass, are listed in Table |l|. In or- 
der to convert the published photometry to a uniform 
set of fluxes we applied our own corrections for the light 
of the probable host galaxy (Holland & Hjorth 1999] ). 
In cases where the light from the galaxy had already 
been subtracted in the published photometry we added 
it back using the magnitude for the galaxy that was given 
in the source then subtracted a uniform estimate of the 
light fro m th e galaxy based on the data of Holland & 
Hjorth ( |1999|) (V = 24.25 ± 0.07, see Table |). Castro- 
Tirado et al. (1999) used UBVRI photometry to determine 



Filter 


A (A) 


/„,o (erg cm -1 s~ x Hz -1 ) 


B 


4448 


4.02 x 10 _aj 


V 


5505 


3.59 x 10~ 20 


Gunn r 


6538 


2.96 x 10" 20 


R 


6588 


3.02 x 10" 20 


I 


8060 


2.38 x 10 -20 



Table 2. The observed magnitudes (from Castro-Tirado 
et al. 1999), and the adopted magnitudes and fluxes for 



the probable host galaxy of GRB 990123. 



Filter 


Castro-Tirado 


adopted 


adopted 




mag 


mag 


flux (/iJy) 


B 


24.23 ±0.10 


24.35 ±0.11 


0.729 ±0.074 


V 


24.20 ± 0.15 


24.25 ±0.07 


0.716 ±0.046 


Gunn r 




24.05 ± 0.07 


0.711 ±0.046 


R 


23.77 ±0.10 


24.07 ± 0.08 


0.709 ± 0.052 


I 


23.65 ±0.16 


23.83 ± 0.08 


0.701 ± 0.052 



a spectral energy distribution of f v oc i/^ sal = constant, 
i.e., /3g a i « 0, for the host galaxy. Therefore, we converted 
the Holland & Hjorth ( |1999| ) V-band flux for the galaxy 
to the other band passes using mag = V + K mag — Ky — 
2.5/3 log 10 (f mag /iv) and (3 ga \ = where Ky and K mag are 
the photometric zero points of the two band passes. The 
Castro-Tirado et al. ( 1999] ) photometry, and our adopted 
(dereddened) magnitudes and fluxes for the probable host 
galaxy, are listed in Table |2| A flat spectrum is consistent 
with the probable host galaxy being a starburst galaxy. 

The Galactic reddening in the direc tion o f GRB 990123 
is E B -v = 0.016±0.020 (Schlegel et al. |l998| ). We assumed 
that there is no extinction in the host galaxy, but this is 
probably not a good assumption since it is unlikely that a 
starburst galaxy is devoid of dust. Reddening due to dust 
in the host galaxy may result in us underestimating the 
flux from the host galaxy, and thus overestimating the flux 
from the OA. To avoid this uncertainty we have fit light 
curves to each band pass separately since the systematic 
error in the flux introduced by the unknown extinction will 
only affect the zero-point of the flux in each bandpass, not 
the slopes of the light curves. 

The fits for each band pass, and each fitting function, 
are given in Table [j| The uncertainties in the parameters 
are the formal la errors in the fit and do not include con- 
tributions from the covariance between the parameters. 
The slopes are not strongly correlated with each other, 
but they are moderately correlated with the location of 
the break. None of the residuals show any time depen- 
dence. 
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Table 3. The parameters of the best-fitting light curves for the photometry of the OA associated with GRB 990123 
in each band pass. Eq. (|l]) is a broken power law, Eq. (|^) is the continuous function of Harrison et al. ( |1999[ ), and 
Eq. (||) is the continuous function of Stanek et al. ( 1999[ ). 





Filter 

r niei 




Oi.2 


tb (days) 




ADOF 


UUr 


/TTK 

w 


D 
ID 


1 11 I r\ i o 

— 1.11 zh U.lz 


— -l.ol ± U.14 


1 C7 _1_ f\ K1 
1.0 ( ± U.Ol 


Ci QO _L A 9C 

y.oz ± 4.oo 


i noo/i 
1.UZZ4 





1 7 
V 


1 OS -4- n A c 
— l.Zo ± U.4o 


— l.Oo ± U.lD 


i 71 l n on 


o no 4-17 


n vori 


o 
o 




Gunn r 


— 1.11 ± 0.12 


—2.01 ± 0.16 


nil I n o o 

2.11 ± 0.83 


6.22 ± 3.94 


2.6654 


4 




R 


-1.17 ±0.30 


-1.57 ±0.11 


1.70 ± 0.22 


10.74 ±3.01 


0.9111 


31 




I 


-1.10 ±0.42 


-1.64 ±0.10 


1.36 ± 0.70 


15.72 ± 10.83 


0.6922 


1 


§ 


B 


-0.83 ± 0.38 


-1.81 ±0.25 


0.76 ± 0.93 


21.89 ±34.83 


0.5445 


6 




V 


-0.76 ± 1.82 


-1.57±0.19 


0.39 ±0.64 


55.90 ± 131.49 


0.8296 


8 




Gunn r 


-0.91 ± 0.24 


-2.11 ±0.31 


1.45 ± 1.34 


9.34 ± 12.57 


3.6172 


4 




R 


-0.73 ± 1.16 


-1.57 ±0.22 


0.54 ±0.84 


40.55 ±71.25 


0.8444 


31 




I 


19.07 ± • • ■ 


-1.61 ±0.06 


0.70 ± 0.80 


26.00 ± 47.72 


0.7975 


1 


i 


B 


-0.75 ±0.70 


-1.90 ±0.43 


0.92 ± 2.11 


17.22 ± 53.45 


0.5647 


6 




V 


-1.12 ± 0.60 


-1.71 ±0.26 


1.40 ±4.29 


11.09 ±48.05 


0.8678 


8 




Gunn r 


-0.87 ±0.38 


-2.25 ±0.47 


1.78 ± 2.46 


6.78 ± 14.37 


3.9109 


4 




R 


-0.56 ±3.23 


-1.63 ±0.68 


0.54 ±2.63 


40.64 ± 214.63 


0.8722 


31 




I 


2.55 ±9.37 


-1.62 ±0.10 


0.54 ±0.24 


36.76 ± 30.45 


0.7451 


1 



The weighted mean values for each parameter, their 
standard errors (SE), and the x 2 P er degree of freedom 
(Xdof) values, are listed in Table for each fitting func- 
tion. The BVrRI data was used to compute the mean 
values of the fits to the broken power-law (Eq. (|l])) and 
the BVrR data was used to compute the mean values of 
the fits to the two continuous functions (Eqs. (||) and (||)). 
The /-band data were not used when computing the lat- 
ter two sets of mean values because the fits to the /-band 
data were not reliable for t < tf,. The mean values of a± 
and «2 agree to within one a regardless of which function 
is fit to the data, but the mean value of tb is significantly 
smaller for Eq. (||) (Harrison's function) than it is for the 
other two fitting functions. However, the broken power 
law (Eq. (Q)) gives considerably less scatter in the val- 
ues of each parameter, particularly tb, for different band 
passes than the other two fitting functions do. Therefore, 
we conclude that the broken power law gives a better fit 
to the GRB 990123 light curves than either of the con- 
tinuous functions do. This suggests that the break in the 
light curve may have been very rapid. The la error in the 
break time for a broken power law is 0.19 days (see Ta- 
ble U), which suggests that the break may have occurred 
over a period of only a few hours. Figs. show the data 
for each band pass with the best-fitting broken power laws 
superimposed. 

To test if the break time depends on the wavelength we 
computed the probabilities that the break times in each 
band pass are consistent with the mean break time. If the 
individual break times are not consistent with the mean 
then this would suggest that the break between the a.\ 
and a2 regimes depends on wavelength, which would be 
inconsistent with the predictions of the relativistic fire- 
ball model. This method makes no assumptions about the 



Table 4. The mean values of ai, ct2, and tb for the best 
fits of Eq. (H) (a broken power law) to the BVrRI data, 
and of Eq. (2|) (Harrison's function) and Eq. (||) (Stanek's 
function) to the BVrR data, for the OA associated with 
GRB 990123. P is the probability that each parameter 
depends on wavelength (see the text). 





Parameter 


Mean ± SE 


2 

X4 


P 


# 


Ql 


-1.12 ± 0.08 


0.0387 


0.0028 




a 2 


-1.69 ±0.06 


1.6616 


0.8442 




t b (days) 


1.68 ±0.19 


0.1331 


0.0297 




Ql 


-0.88 ±0.20 


0.0180 


0.0033 




Q.2 


-1.70 ±0.12 


0.9197 


0.5697 




t b (days) 


0.61 ±0.42 


0.1814 


0.0909 




ai 


-0.91 ±0.29 


0.0658 


0.0220 




Q-2 


-1.83 ±0.19 


0.3748 


0.2288 




t b (days) 


1.11 ± 1.30 


0.0446 


0.0125 



physical origin of the light curve. We checked if the break 
time was dependent on wavelength by evaluating the % 2 
statistic for the weighted mean value of tb- The Xdof 
value, and P, the probability that we can reject the null 
hypothesis that the break time seen in each band pass is 
consistent with the mean break time (i.e. P is the proba- 
bility that the break time changes with wavelength), are 
listed in Table ^. Similar calculations were done for the 
slopes a\ and «2- The P values are all less than 0.9, so we 
conclude that there is no evidence that the break time, or 
the slopes, vary with wavelength. 
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22 

24 

26 

28 
1 

0.5 

3 

2 o 

w 
aj 
K 

-0.5 
1 





1 1 1 1 1 1 




L. 




Xl = 1-0224 2 


|_ a t = - 


-1.11 ± 0.12 ^\ 




: « 2 = - 


-1.81 ± 0.14 ^ 1 


-_ 


f — 1 


j ' i u.iji u.d_y b 




: B (t„) 


= 21.59 ± 0.51 


>» 

1 1 f^H h 


H 1 h 

: RMS = 


— 1 1 1 1 1 1 

0.093 




— 


-Kf-r- 


— -t ^ 




GRB 


990123 B-Band : 




, lii 





-10 12 
log 10 ( Days since 23.4076 Jan. 1999 ) 

Fig. 1. The upper panel shows the best-fitting broken 
power law (Eq. (0)) for the GRB 990123 B-band pho- 
tometry. The horizontal bar shows the ler uncertainty in 
the time of the break. The lower panel shows the residuals 
in the fit. The residuals are defined as (Bat — -Bobs)- The 
uncertainties in the residuals are the uncertainties in the 
observed data. The magnitudes have been corrected for 
Galactic extinction in the direction of GRB 990123, but 
not for extinction in the host galaxy. 

3. 2. The Spectrum of the Optical Transient 

We have assumed that the OA for GRB 990123 has a 
power-law spectrum of the form 



M") 



(4) 



where j v (v) is the flux, in /iJy, at the frequency v, and /?oa 
is the spectral index. The choice of 10 15 Hz as a reference 
frequency is arbitrary, but this value gives v/\Q 15 « 1, 
which results in stable fits of the flux data to Eq. (|4|). 

We estimated the spectrum at a series of times by 
taking the B-, V-, R-, and /-band data in 0.1-day in- 
tervals and fitting it to Eq. (|4|). Only intervals that con- 
tained at least three observations spanning two or more 
filters were considered. Fig. ^| shows the values of /?oa de- 
termined this way. The weighted mean spectral index is 
A^X = -0.750 ± 0.068 (SE) with X j = 0.4711. The low 
Xdof value suggests that we can reject the null hypothe- 
sis, that the spectral index is constant with time, at the 
17% confidence level. This is too low to be able to reject 
the null hypothesis, so we conclude that there is no evi- 
dence for a variable /?oa- Our value of /?oa is con sisten t 
with the direct measurement of Andersen et al. ( 1999|) . 



18 

20 

22 

24 

26 

28 
1 

0.5 




3 
-a 



0) 



-0.5 



X% = 0.7251 



a, = -1.28 ± 0.48 
a 2 = -1.58 ± 0.16 
t b = 1.71 ± 2.20 days 



RMS = 0.128 



GRB 990123 V-Band 



-10 12 
log 10 ( Days since 23.4076 Jan. 1999 ) 

Fig. 2. This figure shows the best-fitting broken power 
law, and the residuals, for the GRB 990123 F-band data. 
The large uncertainty in the time of the break (tb = 1.71 ± 
2.20 days) indicates that the data could be well fit by a 
single power law. 



They found f3 A = -0.69 ± 0.10 between 4000 A < A < 
5700 A from a spectrum of the OA that was taken 0.8 
hours after the burst. 



4. GRB 990510 Data 

4--1- The Photometry 

We collected the published optical observations of 
the OA associated with GRB 990510 from eleve n 
sources: Beuer mann et al. ( 1999 ), Covinoet al. ( 1999 ), 
Fruchter et al. 



Udalski (1999a 1999b, 1999c), and Stanek et al. (1999) 



( |1999aj ), Gal ama et al. (|1999t| ), Harrison 
et al. fll99SD, Marconi et al. fll999a[ |l999b| ), Pietrzyhski & 



These data consists of 182 observations in the Johnson V- 
band, and in the Kron-Cousins R-, and /-bands. The ob- 
servations were taken at seven telescopes up to 40 days af- 
ter the burst. The published magnitudes were converted to 
fluxes in the same manner as described in Sect. 3.1. No cor- 



rections were made for a host galaxy since no host galaxy 
has been observed for this burst to a limiting magnitude of 
R = 27.6 (Beuermann et al. |199S| )F|. 



The Galactic redden- 
ing in the direction of GRB 990510 is E B -v = 0.203±0.020 



1 Recently Bloom (2000) and Fruchter et al (200C) reported 
evidence for a faint (V = 28.5 ± 0.5) host galaxy. This corre- 
sponds to a flux of ~ 0.014 fj,3y in the V band, which is not 
sufficient to affect the results of our fits. 
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Fig. 3. This figure shows the best-fitting broken power 
law, and the residuals, for the GRB 990123 Gunn r-band 
data. 
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Fig. 4. This figure shows the best-fitting broken power 
law, and the residuals, for the GRB 990123 i?-band data. 



(Schlegcl ct al. 1998). We assumed that the extinction in 



the host galaxy was Ay = 0. 

The fits for each band pass, and each fitting function, 
are given in Table ||. The uncertainties in the parameters 
are the formal ler errors in the fits and do not include 
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t b = 1.36 ± 0.70 days 
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Fig. 5. This figure shows the best-fitting broken power 
law, and the residuals, for the GRB 990123 /-band data. 



contributions from the covariance between the parame- 
ters. The residuals for all three fitting functions show some 
structure near the location of the break. This suggests that 
none of the three fitting functions which we used in this 
paper accurately reproduce the temporal behavior of the 
light curve during the transition from the ot\ regime to the 
a-i regime. The last two V^-band observations, and the last 
i?-band observation, were derived from images taken with 
the HST/STIS and were transformed from the 50CCD 
(clear) aperture to the V and R bands assuming a spectral 
-0.6 (Fruchter et al. 



index of /3 A = -0.6 (Fruchter et al. |l999a| ). These ob- 
servations were taken approximately one month after the 
other observations, so the large residuals for these three 
data points may be due to a change in the shape of the 
OA's spectrum during that period, or systematic uncer- 
tainties in transforming the 50CCD magnitudes to VR 
magnitudes. 

The weighted mean values for each parameter, their 
standard errors, and the Xdof values are listed in Table || 
for the fits to each fitting function. The mean values of ax, 
Q.2, and tb agree to within approximately 1.5er regardless 
of which function is fit to the data. However, there are sys- 
tematic variations in the residuals of the fits of Eq. ([!]) (the 
broken power law) to the V- and i?-band data. This indi- 
cates that one of the continuous functions (Eq. (0) or (||)) 
provides a better fit to the data than the broken power 
law (Eq. (|l|)) does, which in turn suggests that the break 
in the light curve was not nearly instantaneous as was 
the case with GRB 990123, but may have occurred over a 
period of approximately one day. Eq. (|^) (Harrison's func- 
tion) gives the smallest residuals, and the least scatter in 
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Table 5. The parameters of the best-fitting light curves for the photometry of the OA associated with GRB 990510 
in each band pass. Eq. (|l]) is a broken power law, Eq. (|^) is the continuous function of Harrison et al. ( |1999[ ), and 
Eq. (||) is the continuous function of Stanek et al. ( 1999[ ). 



Eq. 


Filter 




a 2 


h 


fAh) 


Xdof 


DOF 


© 


V 


-0.78 ±0.12 


-2.08 ±0.20 


1.23 ±0.22 


62.72 ± 13.20 


1.5183 


63 




R 


-0.69 ±0.11 


-1.72 ±0.13 


0.87 ±0.12 


96.06 ± 15.50 


1.2569 


57 




I 


-1.16 ±0.34 


-1.73 ±0.29 


1.33 ±0.66 


58.67 ±41.41 


0.9355 


50 


© 


V 


-0.62 ±0.19 


-2.16 ±0.29 


1.05 ±0.47 


65.54 ± 37.73 


0.6395 


63 




R 


-0.46 ±0.20 


-1.85 ±0.26 


0.70 ± 0.35 


114.54 ±63.24 


0.7292 


57 




I 


-0.67 ± 1.79 


-1.84 ± 1.22 


0.80 ± 2.26 


106.44 ±351.17 


0.8679 


50 


© 


V 


-0.58 ± 0.23 


-2.32 ±0.37 


1.25 ±0.71 


50.70 ± 39.84 


0.4218 


63 




R 


-0.40 ±0.29 


-2.00 ±0.36 


0.81 ±0.60 


96.02 ± 84.48 


0.9763 


57 




I 


-0.51 ±2.08 


-2.02 ± 1.16 


0.98 ± 2.85 


83.55 ± 310.43 


0.8547 


50 



0.5 



-0.5 



■1.5 



R = -0.750 ± 0.068 



GRB 990123 



12 3 

Time (Days since 23.4076 Jan. 1999) 

Fig. 6. This figure shows /3oa, the spectral index of the 
light from the OA associated with GRB 990123 (filled cir- 
cles), as a function of time. The error bars are defined 



as a/3 



12iLi ( a f,i/ ' fi) 2 where aj^ is the uncertainty 



in each flux value, fi. The solid line shows the weighted 
mean value of /3oa while the dashed lines show the la 
uncertainty in /?oa- The open square shows the Andersen 
et al. ( |1999D spectral value. The data are consistent with 
a constant /3oa of —0.750 ± 0.068. There is no evidence 
that /3qa varies with time. 



the values of the parameters in different band passes, so 
we have adopted it as the best fitting function. The P val- 
ues in Table |^ are all less than 0.7, so we conclude that 
there is no evidence for any variation in u\ , cti , or tb with 
wavelength. Figs. show the data for each band pass 
with the best-fitting continuous functions superimposed. 



Table 6. The mean values of ai, ct2, and tb for the best 
fits of Eq. (|]J) (a broken power law), Eq. (|^) (Harrison's 
function), and Eq. (|3|) (Stanek's function) to the VRI data 
for the OA associated with GRB 990510. P is the prob- 
ability that each parameter depends on wavelength (see 
the text). 





Parameter 


Mean ± SE 


xi 


P 


# 


Ql 


-0.75 ±0.18 


1.8114 


0.5957 




a 2 


-1.81 ±0.10 


1.1878 


0.6951 




tt 


0.96 ±0.10 


1.1913 


0.6962 


(!) 


ai 


-0.54 ±0.14 


0.1707 


0.1569 




OL1 


-1.98 ±0.19 


0.3239 


0.2767 




tb 


0.82 ±0.28 


0.1784 


0.1634 


® 


a?i 


-0.51 ±0.18 


0.1182 


0.1115 




a L > 


-2.15 ±0.25 


0.1986 


0.1801 




U 


0.99 ± 0.45 


0.1120 


0.1060 



4-.2. The Spectrum of the Optical Transient 

We estimated the spectral index for the OA associated 
with GRB 990510 in the same way as we did for the OA 
associated with GRB 990123 (see Sect. |J). Fi S- " 



|Hj shows 

the values of /?oa that were determined in this manner. 
The weighted mean spectral index is (3oa = — 0.531 ± 
0.019 (SE) with X14 = 1-6598, which suggests that we can 
reject the hypothesis that the spectral index is constant 
with time at the 95% confidence level. Our mean value for 
/3oa is consistent (within 1.5a) with the value obtained by 
Beuermann et al. ( 1999] ) using long-slit spectrophotometry 
3.9 days after the burst. They found [3oa = — 0.55±0.10 
in the interval 4900 A < A < 9000 A. 

Our result is consistent with the spectral index increas- 
ing from /3oa ~ at the time of the burst to /?oa ~ —0.6 
at the time of the break in the light curve (tb = 0.82 ±0.28 
days after the burst). After the break the spectral index 
stays constant at /3oa of —0.62 ± 0.06 , which is consis- 
tent with the Beuermann et al. ( |l999| ) result. We find 
(3qa = —1.29 ±0.23 3.6 days after the burst. This is incon- 
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Fig. 7. The upper panel shows the best-fitting continuous 
function (Eq. (§), Harrison's function) for the GRB 990510 
V^-band photometry. The horizontal bar shows the la un- 
certainty in the time of the break. The lower panel shows 
the residuals in the fit in the sense (Vet — Vobs)- The mag- 
nitudes have been corrected for Galactic extinction in the 
direction of GRB 990510. 
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Fig. 8. This figure shows the best-fitting continuous func- 
tion (Eq. (|)), and the residuals, for the GRB 990510 R- 
band data. 



Fig. 9. This figure shows the best-fitting continuous func- 
tion (Eq. (§)), and the residuals, for the GRB 990510 I- 
band data. The large uncertainty in the time of the break 
(if, = 0.80 ± 2.26 days) indicates that the data can be well 
fit by a single power law. 



sistcnt with the Bcucrmann et al. (1999) data, but is con- 
sistent with the spectral index continuing to increase after 
the break in the light curve occurred. However, we believe 
that the direct measurement of Beuermann et al. (1999) 
is more reliable than our indirect measurements of /?oa 
so we conclude that it is unlikely that the spectral index 
continues to increase after the break. 



5. Models 

The generic model for an OA associated with a GRB is a 
fireball expanding relativistically into an ambient medium 
of number density n, and decelerating as it sweeps up mat- 
ter. The shock between the fireball and the medium accel- 
erates electrons to relativistic energies, and gives them a 
power-law distribution of energies, N(~f) oc 7 _p . The total 
energy in the electrons is parameterized by e e , the ratio 
between the energy in electrons and the energy in nucle- 
oids. Similarly, the magnetic field strength is parameter- 
ized by €b , the ratio of the comoving field energy density 
and the post-shock nucleon energy density (e.g. Wijers & 
Galama|i99S|). 



At a given instant in the fireball evolution the instan- 
taneous spectrum of the OA is determined by the local 
conditions in the shock region. In particular, the electron 
energy index p determines the spectral slope of the syn- 
chrotron emission f v oc t a v^ , e.g., through the relation 
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For a collimated outflow in an ambient medium of con- 




1 2 3 

Time (Days since 10.432 May 1999) 

Fig. 10. This figure shows /3oa, the spectral index of the 
light from the OA associated with GRB 990510 (filled 
circles), as a function of time. The solid line shows the 
weighted mean value of /3oa while the dashed lines show 
the la uncertainty in /3oA- The open square shows the 
Beuermann et al. ( 1999) spectrophotometric data. The 



data are inconsistent with a constant /3oa at the 95% 
confidence level (92% confidence level if the spectropho- 
tometric datum is included). 



/3 = -(p- l)/2. Following Sari et al. (119981) , we have, for 
the fast cooling regime, 



fu 



{v/v c )- 1/2 f Vml 



V c < V < Vn 



(v m /v c ) 1/2 {v/v m ) p/2 f Vml V>V C 
while in the slow cooling regime we have 



(5) 



{v/v v 



-(p-l)/2 



< V < v c 



{VolVmT {p - l), \v/VcT P,2 f» m , Vc<V- 



(6) 



The synchrotron peak frequency, v my the flux at that fre- 
quency f Vm , and the frequency v c corresponding to the 
electron Lorentz factor that separates the fast cooling elec- 
trons from the slow cooling ones, all depend on the model 
parameters and the fireball geometry (spherical or colli- 
mated). Explicit expression for these cases can be found 
in Sari et al. (|1998| ), Wijers & Galama ( |1998| ) (spherical), 
and Rhoads (1999) (jet). We do not consider the spectrum 
below v c and v m in the fast and slow cooling regimes, re- 
spectively, as this is not of interest in the present discus- 
sion. The time dependence of v m , u c , and f Um is different 
in spherical and collim ated models and h as bee n discussed 
by e.g. Rhoads ( 1999 ; see also Sari et al. 1999 ). 



stant number density n, we have (Rhoads 1999) 



3.66 ■ I0 16 x v eie 



2 1/2 



E, 



r>2 



1/2 



1 + 2 



,V2 



.3/2 



Hz. 



1.04- 10 n e 



-3/2^-2/3 - 
3 E 52 n 



5/6 Hz, and 



f Vm = 1.8 • 10" 



(1 + z)ej <t> P y-gr 



,1/2. 



4.82 Gpc 



(7) 



(8) 



Jy.(9) 



Here, x p is the dimensionless frequency at maximum flux 
and 4> p is the dimensionless flux at that frequency (Wijers 
& Galama 1998). These are determined by the index p 
of the electron energy distribution. The burst energy in 
units of 10 52 erg is denoted by E 52 , is the half angle 
of the collimated jet, td is the time in days, and is the 
luminosity distance to the burst. 

To connect the optical observations in a given band 
pass directly to the model parameters we rewrite, for ex- 
ample, Eq. (t|) as 



(10) 



for v m < v < v c in the slow cooling regime, while for 
v > v c we would have 



y p/2 U 



(p-l)/2 



,p/2 



(11) 



Similar expression may be obtained in the fast cooling 
regime. Combining v m and f„ m according to Eq. (|l0|), or 
i> m , v c , and f Vm according to Eq. (pi]), depending on the 
frequency range of interest, allows us to directly connect 
the observations at a given frequency to the model pa- 
rameters. The general formula that is obtained this way 
is complicated but may be simplified in specific cases (see 
Eq. 

The properties of the observed light curve of the OA, 
such as the rate of decay a, are primarily determined by 
the hydrodynamic evolution of the fireball. In particular, 
spherical and collimated bursts give rise to different forms 
of light curves. The adiabatic evolution of the Lore ntz fac- 
tor, r , of th e fireball is given by (Sari et al. 1998 ; Sari & 
Piran |l999| ) 



T = 6.65 



\ V8 -3/8 

U ) d ' 



(12) 



A spherical outflow generally results in a uniformly de- 
caying light curve that can have modest breaks which are 
due to the different time evolution of the characteristic fre- 



quencies v m and v c (Sari et al. 1998). A pronounced break 



in the light curve may be observed in a spherical outflow 
that becomes non-relativistic, but requires the ambient 
medium to be denser, by several orders of magnitude, than 
generally expected (Wijers et al. 1997 ; Dai & Lu 1999 ). 
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The steepening of the light curve in the case of a spherical 
fireball entering a non-relativistic phase can be shown to 
be Aa = (a± — a%) = — {a\ + 3/5). 

Collimated outflows, on the other hand, give rise to a 
light curve that resembles a broken power-law which be- 
comes steeper between several hours and several days after 
the gamma-ray event. There are two general possibilities. 

1. There is a collimated outflow with a fixed opening an- 
gle 8 n . When T decreases to T w l/#o> the light curve 
steepens. This steepening occurs because after the out- 
flow has decelerated the observer only receives radi- 
ation e mitted within the colli mated beam (Meszaros 
k Rees |200C|; Kulkarni et al. |l999a|). The amount of 



steepening is Aa = 3/4. This is a purely geometrical 
effect, so the intrinsic properties of the fireball are not 
affected by the transition to the T < 1/8q regime. 
The collimated outflow may expan d sideways as 
well as radially (e.g. Rhoads 1999; Panaitescu & 



Meszaros 1998). Initially, the relativistic beaming dom- 
inates over the sideways expansion, but when side- 
ways expansion sets in, the radial expansion essentially 
comes to a halt, which results in a steepening of the 
light curve. Simple analytic estimates show that this 
steepening is Aa = 1 — ai/3 (Rhoads 1999| ), while de- 
tailed numerical simulations show that the steepening 
is expected to be more gradual both in time and slope 
(Moderski et al. |2000p. 



6. Interpretation 

6.1. GRB 990123 

The synchrotron peak frequency, v m , was most likely 
in the radio already on 1999 January 24.65 (Galama 
et al. 1999a ), implying that the electrons were already 
in the slow cooling regime. The flux at frequency v would 
then be given by Eq. (|Io|). The isotropic energy of the 
burst, £52, that was estimated from the gamma-ray data, 
is -E52 = 338 ± 1 at the observed redshift of z = 1.6. 

The spectral index that we determined for 
GRB 990123, f3 = -0.750 ± 0.068, corresponds to 
an electron index of p = 2.500 ± 0.136. There is no 
evidence for variations in f3 across the break in the light 
curve, so we can conclude that p is constant during the 
observed evolution of the light curve. The reliability of 
P and p obtained in this way is questionable, however, 
as the effect of extinction on (3 is unknown. The light 
curve decay indices, ai and a2, on the other hand, are 
also a function of p, and are not affected by extinction. 
Therefore, we have elected to use our fits to the measured 
light curves to obtain the electron index, and use that to 
infer the spectral index. 

In Table | we compare the electron energy index that 
we calculate from the light curve indices using three differ- 
ent models and assuming that the ambient medium has a 
constant density. During the initial phase, all models look 



spherical to the observer, hence p is the same for all of 
them. After the break in the light curve, different models 
give different relationship between the light curve decay 
rate and p, as indicated in Table |7|. The last line in Table [?] 
gives the size of the break that is predicted by each model 
and should be compared to the observed size of the break, 
Aa = 0.57±0.10. A spherically expanding fireball entering 
a non-relativistic regime at the time of light curve break 
gives the best o verall agreem ent w ith the data (see also 
Dai & Lu |1999| ; Wang et al. |200C|) . That interpretation, 
however, does have problems in accounting simultaneously 
for t he optic al and radio properties of the burst (Kulkarni 
et al. 1999b ). It is worth noting that the sideways expand- 
ing jet predicts a break in the light curve of Aa = 1.37, 
which is more than twice the size of the observed break. 

Our preferred model in interpreting this burst is that 
of a collimated outflow of a fixed opening angle. The 
main reason for this choice is the model's ability to self- 
consistently account for the observations at all wave- 
lengths. If we accept the constant #0 model as the best 
description of GRB 990123, we can deduce the remaining 
model parameters. Evaluating Eq. (jl|) at td = h = 1.68, 
and using E52 = 338, we find that the Lorentz factor of 
the fireball at the time of the break is 



T = 5.47 



',2 



lln- 1 ' 8 . 



(13) 



Therefore, interpreting the break in the light curve as a 
geometrical effect, we find the opening angle of the colli- 
mated outflow to be 



1 
f 



1 

IT' 



*V8 w 5V/ 8 , 



(14) 



which enables us to eliminate the opening angle from the 
expressions for v rn and /„ m (Eqs. (0) and @). The open- 
ing angle is weakly dependent on the number density of 
the medium where the burst occurred. Typical number 
densities in star-forming regions, for example, can vary 
from n w 1 to ji ^ 1000. However, this range of number 
densities corresponds to opening angles of 5° 8 ^ 12°, 
so the amount of collimation does not depend strongly on 
the density of the ambient medium. 

The fireball evolution is independent of colour, so ob- 
servations in different bands passes should give identical 
results when evaluating the left hand side of Eq. (|To|) . Us- 
ing the data for Eq. ([!]) (the broken power law) in Table || 
at the time of the break in the corresponding light curve, 
averaging over all five band passes, and rearranging, we 
obtain 



,p-l,(p+l)/4 p (p+l)/4 ( p _9)/32 
fc e t B -°52 n 



1.34-10" 



(15) 



We have retained p in the exponents of the model param- 
eters to show their explicit appearance, but used p — 2.37 
(the mean of p\ and pi for the constant 9q model), to 
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Table 7. The electron energy index before, p\, and after, p2, the break in the light curve, and the size of the break, 
Aa, inferred from the observations of GRB 990123 for the three models discussed in the text. The observed values are 
given for comparison. The observed value of p is computed from our observed value of [3oa = —0.750 ± 0.068. 





Spherical ± non- 


relativistic 


8o constant 


Sideways expansion 


Observed 


Pi 


1-4qi/3 = 


2.49 ±0.11 


l-4ai/3 = 2.49 ±0.11 


l-4ai/3 = 2.49 ±0.11 


2.50 ±0.14 


P2 


(21 - 10q 2 )/15 = 


2.53 ±0.04 


-4a 2 /3 = 2.25 ±0.06 


-a 2 = 1.69 ±0.06 


2.50 ±0.14 


Aa 


-(ai+3/5) = 


0.52 ±0.06 


3/4 


l-ai/3 = 1.37 ±0.02 


0.57 ±0.10 



evaluate the right hand side. Note that the above equa- 
tion constrains four independent model parameters, but 
the isotropic energy of the burst, -E52 = 338 ± 1, can be 
determined directly from observations, so only three pa- 
rameters remain undetermined. 

An additional expression can be obtained using the 2 — 
10 ke V X-r ay data obtained on 1999 January 24.2(Heise 
et al. 1999 ). At that time the cooling frequency was less 
than the frequency of the X-rays [y > u c ), so Eq. ( |il"l) 
applies. Evaluating the left-hand side of Eq. ([ll]) using 
the A-ray observations we find 



£r l e (p-5)/4 ri _ b _9)/32-5/6 jE ;(3p-5)/12 = ^ 



and combining the optical and A-ray results gives 



e 3 R /2 n 5/6 ps 1.33 • 10~ 7 . 



(16) 



(17) 



This gives e B ps 2.6 • 10 -5 and e e ps 0.79 if n = 1 cm" 3 , 
or e B ps 5.6 • 10~ 7 and e e ps 2.94 if n = 1000 cm" 3 . An 
additional relation is needed in order to determine n in- 
dependently. 

An alternative, but related interpretation is collimated 
outflow in an ambient medium with a number density 
distribution of the form n(r) cx r~ s (see Panaitescu 
et al. |1998| , Meszaros et al. |1998[ ). Here, 6 = repre- 
sents a homogeneous distribution, while 6 = 2 repre- 
sents a distribution characteristic of a pre-existing stel- 
lar wind. The light curve decay rates will be steeper in a 
wind distribution than if the density is constant, and the 
break in the light curve will be less pronounced. In general 
Aa = (3 — 6)/{A — 6). This reproduces Aa = 3/4 when 

5 = 0, as discussed above, but gives Aa = 0.5 for a stellar 
wind density distribution (6 = 2). 

Interpreting the burst using this model we obtain 

6 = 1.67 ± 0.54, which is consistent with a stellar wind 
density distribution. This model adds 6 as a new parame- 
ter, but constrains /?, since the spectral index is now fixed 
by the light curve decay rates (a\ and 02) and 6. We 
obtain = —0.51 ± 0.14, which is marginally consistent 
with the spectral index inferred from the data, resulting 
inp = 2.02 ±0.28. The other model parameters, discussed 
above are similar to what is found in the 6 = case. 



6.2. GRB 990510 

GRB 990510 was a strong burst, although none of its 
observed properties were extraordinary in any way (Wi- 
jers et al 1999| ). The lower limit to its rcdshift is z = 



1.619 ±0.002 (Vreeswijk et al. |1999| ), and the isotropic en- 
ergy release was -E52 = 17 ± 1. The afterglow was the first 
to show a smooth break in the light curve rather than a 
sharp break at a fixed point in observer time. 

The measured light curve decay index before the break, 
ai = —0.54 ± 0.14, is too small to be explained by adia- 
batic fireball scenarios as it leads to a value for the electron 
energy index, p < 2, that is too small. It is, however, con- 
sistent with a radiative fireball (e.g. Sari et al. [1998 ). In 
that case the light curve is predicted to decay as i -4 / 7 
with a corresponding spectral index of (3 = —1/2. These 
values are both independent of p, and are in good agree- 
ment with the observed value of a±, and the value of f3 
that we inferred from the light curves. After the initial 
radiative stage, the fireball evolves to the more commonly 
observed adiabatic phase. 

In Table |^ we compare the electron energy index, p, 
calculated from the observed value of a\ for the three 
models considered in the previous section. As already 
mentioned the initial slow decay is incompatible with 
the requirement that the electron population have a fi- 
nite energy, i.e., electron energy indices of p < 2 corre- 
spond to infinite energy, and only the sideways expand- 
ing jet model is consistent with the behavior of the light 
curve at late time. This model also gives good agree- 
ment for the size of the break in the observed light curve, 
Aa = 1.44±0.17. The smooth break in the light curve (see 
also Stanek et al. 199G, Israel et al. ( 1999| ), and Harrison 
et al. 1999), further supports that interpretation (Mod- 
erski et al. 200C). The late time spectral index, calcu- 
lated from P2, is /3oa = —0.49 ± 0.10. This is in good 
agreement with the value determine from spectrophotom- 



etry (Poa = — 0.55 ± 0.10, Beuermann et al. 1999), and 
the mean value that we calculated from the light curve 
(Poa = —0.53 ±0.02). Therefore, the sideways expanding 
jet, with a radiative initial phase, is our preferred model. 

Applying an analysis similar to the one that we did for 
GRB 990123, we obtain an opening angle of 9q ps ^n 1 ^ , 
which is con sisten t with the opening angle derived by Har- 
rison et al. ( 1999 ) for GRB 990510 and the opening angle 
that we derived for GRB 990123. In principle it is possible 
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to derive e e and es in the same manner that we did for 
GRB 990123. However, to do this we need an observation 
of the flux at v > v c , and, unfortunately, there are no ob- 
servations in the literature at these frequencies. Therefore, 
we are not able to break the degeneracy between e e , e^, 
and n that is seen in Eq. (p"5|). 

7. Conclusions 

It is possible to learn a great deal about the nature of 
a GRB, and the environment where the burst occurred, 
from the optical light curve of the burst's OA. However, it 
is important that the decaying light curve be well-sampled 
in order to accurately determine the time and shape of the 
break, as well as the slope of the light curve before and 
after the break. 

GRB 990123 and GRB 990510 are the first two GRBs 
for which breaks have been observed in their optical light 
curves. We have collected all of the published photome- 
try for these two bursts and determined the time of the 
break and the slopes of these light curves before and af- 
ter the break using three different fitting functions. We 
also estimated the spectral indices of each OA from near- 
simultaneous broad-band photometry. These parameters 
were used to constrain models for the nature of each burst. 
Our results suggest the following. 

1. Both GRB 990123 and GRB 990510 exhibit breaks 
in their optical light curves at times between approx- 
imately one and two days after the bursts. There is 
no evidence that the time of the break depends on fre- 
quency. The break in GRB 990123's optical light curve 
appears to have occurred very fast, perhaps over a pe- 
riod of less than five hours. This is in sharp contrast to 
the break in the optical light curve for GRB 990510, 
which appears to have occurred slowly over approxi- 
mately one day. 

2. Unlike GRB 000301C (Jensen et al. |2000D the parame- 
ters of the light curve «2, and i&) for GRB 990123 
and GRB 990510 are not strongly dependent on the fit- 
ting function. This means that our interpretations are 
not strongly dependent on the type of fitting function 
used to derive the slopes and break time. However, the 
choice of fitting function can influence the interpreta- 
tion of the light curve for some GRBs. 

3. Our favoured interpretation for the optical light curve 
of GRB 990123 is that the afterglow was collimated 
with a fixed opening angle of 9 « 5°. Therefore, the 
observed break in the light curve is a geometric effect 
that occurs when the relativistic expansion of the fire- 
ball slows to r « l/#o • This model reduces the energy 
released by the GRB by a factor of » 300 to E52 =1.1, 
which corresponds to the conversion of rs 0.01 Solar 
masses of material to gamma radiation. The magnetic 
field strength is es ~ 2.6 x 10~ 7 if the local number 
density is k 1 cm -3 and cb ~ 5.6 x 10 -5 if the local 



number density is ~ 1000 cm 3 . This is well below the 
equipartition energy. 
4. The OA's light curve for GRB 990510 is consistent 
with a collimated outflow that is initially radiative. 
The transition to an adiabatic phase, and the onset of 
a sideways expansion of the jet, occurred near the time 
of the break. The opening angle was 9q ks 5° at that 
time, the same as GRB 990123. We are unable to es- 
timate reliably the magnetic field strength because no 
measurements are available of the flux at frequencies 
above v c , the electron cooling frequency. 

We have investigated several models for GRB 990123 
and GRB 990510 and find that the OAs for these two 
bursts can not be explained with a single model. This 
suggests that the local environment where a GRB occurs 
plays an important role in determining the evolution of 
the fireball, and the observed light curve of the OA. We 
note, however, that the uncertainties in determining the 
break time and slopes of the light curves make it difficult 
to unambiguously determine the physics of the optical af- 
terglow and its environment. 
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